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ABSTRACT 

In this paper we are going to discuss and analyze the different methods which are developed to detect the 
Lung nodules which cause the lung cancer. At the end of analyzing different methods, the new methodology 
of detecting the lung nodules using Spline Wavelet technique has been proposed in this paper. Continuous 
modeling of data often required in medical imaging, Polynomial Splines are especially useful to consider 
image data as continuum rather than discrete array of pixels. The multi resolution property of Splines 
makes them prime candidates for constructing wavelet bases. Wavelet tool also let us to compress the origi- 
nal CT image to greater factor without any sacrifice in accuracy of nodule detection. Different Algorithms 
for segmentation/ detection of lung nodules from CT image is discussed in this paper. (Int J Biomed Sci 2013; 
9 (1): 9-17) 
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INTRODUCTION 

Lung cancer is one of the serious cancer which caus- 
ing more death world wide. Statistics says 28 percent 
of overall cancer deaths are due to lung cancer. It is the 
number one cause of death from cancer every year and 
the second most diagnosed after breast and prostate can- 
cers (for women and men, respectively). Lung cancer is 
usually found in older persons because it develops over a 
long period of time. Lung cancer is usually visible as small 
round lesion called 'nodules' through Medical images like 
CT. However the major problem in identifying the nodules 
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with CT image is that the characteristics of nodules are 
similar to the characteristics of blood vessels and bronchi 
in terms of size, shape and density. Hence it is necessary 
to enhance and segment the nodules using some special 
Image processing techniques. Early detection of lung can- 
cer using the CT image has the largest chance of saving 
the patient's life. It is important to enhance and detect the 
nodules in CT images in order to identify the lung cancer 
at early stage. Different methods and algorithms are devel- 
oped to effectively detect the nodules. 

In this paper the various methods are discussed and 
one new method is proposed for efficient nodules detec- 
tion. 

PREVIOUS RELATED WORKS 

Initially some works had been carried out to detecting 
the nodules with the help of X-ray (1). X-ray may show a 
mass in the lungs or enlarged lymph nodes. Sometimes 
the chest x-ray is normal, and further tests are needed look 
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for a suspected lung cancer. Even if a mass is found, these 
are not always cancerous and further studies are needed. 
Hence researchers start to work with CT scan image which 
have more information about the tissues of body part com- 
pare to X-ray Some other medical images like MRI and 
PET may have much more useful information than CT, but 
because of the cost factor CT lead the diagnosis process. 
Hence whatever advancement and enhancement carried 
on CT images will help the people greatly. In this section 
we study some of the research carried out with the CT im- 
ages of lungs to find the lung nodules. 

3-D Vascular method of Lung Nodule Classification 

This method includes the enhancement of lung struc- 
tures followed by series of segmentation methods to ex- 
tract the nodule and to form a 3D configuration at an area 
of interest (nodule) (2). The vascular index, aspect ratio, 
circularity, irregularity, extent, compactness and convex- 
ity were also computed as shape features for quantifying 
the nodule boundary. The main disadvantage of the 3D 
methods are it takes more computing time. 

Enhancement and Detection of lung modules with 
Multi-scale filters 

In this method Takemura, Han, chen, Ito and Nishikwa 
(3) proposed a two step method for lung nodules detection 
in CT images. Firstly they used 2D multi-scale filter to de- 
tect the candidates of lung nodules on the slice images and 
then reduce the most of false positive nodules using logical 
AND operator of continuous slices (Fig. 1). 

They take a slice of an original CT image, which has a 
lung nodule, as shown in Fig. 2a, therein, red circle shows 
the location of the true-positive nodule. In Fig. 2a, the en- 
hanced nodules by a 3D multi-scale filter, a 2D multi-scale 
filter and the proposed method are shown in Fig. 2b-2d, 
respectively. 



CT Images 




2D Multi-scale 
filter 




Condidates of 
Lung nodules 
















Reducing most of 
False-positive candidates 




Logical AND operator 
of continuous slices 



In Fig. 3, they showed another example. From Fig. 2 
and Fig. 3, it is obvious that the results of 2D multi-scale 
filter retain too much false positive nodules. The results of 
3D multi-scale filter and this proposed approach got com- 
parable good results. However, 3D multi-scale filter takes 




Figure 2. a, Original CT Image; b, Enhanced nodules by 3D 
multi-scale filter; c, 2D Multi-scale filter; d, Takemura's pro- 
posed method. 




USE 



(c) 



(d) 



Figure 1. Two step method for lung nodules detection in CT 
images 



Figure 3. a, Original CT Image; b, Enhanced nodules by 3D 
multi-scale filter; c, 2D Multi-scale filter; d, Takemura's pro- 
posed method. 
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so much time and it is not applicable in real system, while 
this method can only need the time similar to 2D multi- 
scale filter method. In case of 31 slices, window size is 21 
and scale 4, 3D multi-scale filter takes 1.2 day to process, 
while this proposed method takes only 15 min. 

The objective of this method proposed by Anitha and 
Sridhar (4) is to develop a segmentation system in order 
to assist the surgeons to remove the portion of lung for 
the treatment of certain illness such as lung cancer, and 
tumors. 

This method having following steps: 

a. Image Enhancement and Segmentation. 

b. Fissure Detection for identifying the fissure regions. 

c. Fissure Refinement using region growing method. 

d. Nodule Segmentation using Adaptive Threshold 
The system architecture is presented in Figure 4. 

In this method initially pre-processing of a CT image 
is done to remove the noises present in it. Then the vessels 
are enhanced using the morphological operators followed 
by lung segmentation. The fissure regions are identified, 
enhanced and verified. Lung nodules are segmented using 
adaptive threshold. Through the project they have devel- 
oped a segmentation algorithm for identifying the fissures 
and nodules from CT images. Using the statistical analysis 
such as mean difference between the automatic and man- 
ual segmentation, the results are verified. These statistical 
results are tabulated in Table 1. 

SPLINE WAVELETS 

Researchers are now faced with an ever increasing va- 
riety of wavelet bases to choose from. While the choice of 
the "best" wavelets obviously application-depend, it can 



be useful to isolate a number of properties and features 
that are of general interest to the user. The purpose of this 
section is to present a list of arguments in favor of splines, 
which are unique in a number of ways. 

Wavelets can be classified in four categories: orthogo- 
nal (Battle-Lemarie), semi-orthogonal (e.g., B -spline), 
shift-orthogonal, and biorthogonal (Cohen-Daubechies- 
Feauveau). Unlike most other wavelet bases, splines have 
explicit formulae in both the time and frequency domain, 
which greatly facilitates their manipulation. They allow 
for a progressive transition between the two extreme cases 
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Figure 4. system architecture. 



Table 1. Results for identifying the fissures and nodules from CT images 




Manual 


Automatic 
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of a multiresolution (5): Haar's piecewise constant repre- 
sentation (spline of degree zero) versus Shannon's band- 
limited model (which corresponds to a spline of infinite 
order). Spline wavelets are extremely regular and usually 
symmetric or anti-symmetric. They can be designed to 
have compact support and to achieve optimal time -fre- 
quency localization (B-spline wavelets). The underlying 
scaling functions are the B-splines, which are the short- 
est and most regular scaling functions of order L. Finally, 
splines have the best approximation properties among all 
known wavelets of a given order L. In other words, they 
are the best for approximating smooth functions. 

Spline wavelet in Medical Image processing 

Polynomial splines are especially useful when one 
wishes to consider image data as a continuum rather than 
a discrete array of pixels. Such a continuous modeling of 
the data is often required in medical imaging. Interpola- 
tion, in particular, plays a crucial role at various stages 
of processing. For instance, it is present-explicitly or not- 
for tomographic reconstruction, irrespective of the type 
of algorithm used (filtered back-projection, inverse Fou- 
rier or iterative reconstruction). Another important area 
is medical image visualization; this involves simple 2D 
operations such as image zooming, panning, rotation, or 
3D manipulations, such as reslicing or maximum intensity 
projection 29, which are often used by radiologists. Inter- 
polation models are also required for performing various 
types of image registrations 10, 47; these include intra- 
modal registration for rigid-body motion compensation, 
inter-modal registration of CT, PET and MR data sets of a 
same subject. Considering an image as a continuously-de- 
fined function is also often desirable for feature extraction, 
in particular, contour detection. These are all examples of 
medical imaging tasks that can benefit from the use of 
splines; a more complete inventory is given in Table 2. 

Our purpose is to justify the use of splines in imaging 
applications, emphasizing their ease of use, as well as their 
fundamental properties. Modeling images with splines is 
painless: it essentially amounts to replacing the pixels by 
B-spline basis functions, which are piecewise polynomials 
with a maximum order of differentiability. 

The spline representation is flexible and provides the 
best cost/quality tradeoff among all interpolation meth- 
ods: by increasing the degree, one shifts from a simple 
piecewise linear representation to a higher order one that 
gets closer and closer to being bandlimited. Multiresolu- 
tion properties of splines make them especially attrac- 
tive for multi-scale processing. On the more fundamental 



front, we will show that splines are intimately linked to 
differentials; in fact, the B-splines are the exact mathemat- 
ical translators between the discrete and continuous ver- 
sions of the operator. This is probably the reason why these 
functions play such a fundamental role in wavelet theory. 
Splines may also be justified on variational or statistical 
grounds; in particular, they can be shown to be optimal for 
the representation of fractal-like signals. 

Mathematical model of Polynomial spline 

A polynomial spline (6) of degree n is made up of poly- 
nomial segments of degree n that are connected in a way 
that garantees the continuity of the function and of its de- 
rivative up to order n-l. The joining points between the 
polynomial segments are called knots. In the context of 
the wavelet transform, the knots are equally-spaced and 
typically positioned at the integers. One can thus define a 
hierarchy of spline subspaces of degree n, {Vi n } ieZ, 
where VP is the subspace of L2-functions that are (n-Y) 
times continuously differentiable and are polynomials of 
degree n in each interval (2'k, 2'(k+l)), keZ. The spacing 
between the knot points 2' is controlled by the scale index 
i. Clearly, a function f(x)eViO n that is piecewise polyno- 
mial on each segment (2 10 k, 2 10 (A:+1)) is also included in 
any of the finer subspaces FiPwith ieiO. Thus, we have the 
following inclusion property 

{0}C---F 1 n C---F 0 n C F;"c-cX 2 

Furthermore, it is well known that one can approximate 
any L2-function by a spline as closely as one wishes by let- 
ting the knot spacing (or scale) go to zero ( - infinity). 
This means that the above sequence of nested subspaces 
is dense in L2 and therefore meets all the requirements for 
a multiresolution analysis of LI in the sense defined by 
Mallat. This implies that it is indeed possible to construct 
wavelet bases that are polynomial splines. The best way to 
proceed is to use Schoenberg's representation of splines in 
terms of the B-spline basis functions. In order to satisfy 
the multiresolution inclusion property for any degree n, we 
will use the so-called causal B splines which can be con- 
structed from the («+l)-fold convolution of the indicator 
function in the unit interval (causal B-spline of degree 0) 

n f x 0 0 , s 

cp (x) = cp * . . . * cp (x) 

(n+1) times 

where 

°(x) = ( 1 ' °" X - 1 
\ 0, otherwise 
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Table 2. Summary of medical image application using spline wavelets 



Image processing task 

Tomographic reconstruction 



Sampling grid convertion 



Visualization 



Geometrical correction 



Registration 



Specific operation 

• Filtered backprojection 

• Fourier reconstruction 

• Iterative techniques 

■ 3D + time 

■ Polar-to-cartesian coordinates 

■ Spiral sampling 

• k-space sampling 

■ Scan conversion 
2D operations 

■ Zooming, parnning, rotation 
•Re-sizing, scaling 

■ Stereo imaging 

■ Range, topography 
3D operations 

■ Re-slicing 

■ Max. intensity projection 

• Simulated X-ray projection 
Surface/volume rendering 

■ Iso-surface ray tracing 

• Gradient-based shading 

■ Stereogram 

■ Wide-angle lenses 

• Projective mapping 
■Aspect ratio, tilt 

• Magnetic field distortions 

• Motion compensation 

• Image subtraction 

• Mosaicking 

• Correlation-averaging 

• Patient positioning 

• Retrospective comparisons 

• Multi-modality imaging 

• Stereotactic normalization 

• Brain warping 

■ Contours 

• Ridges 

• Differential geometry 
Contour extraction 

• Snakes and active contours 



Imaging modality 

Commercial CT (X-rays) 
EM 

PET, SPECT 

Dynamic CT, SPECT, PET 
Ultrasound (endovascular) 
Spiral CT, MRI 
MRI 



All 

Fundus camera 
OCT 

CT, MRI, MRA 



CT 
MRI 

Endoscopy 
C-Arm fluoroscopy 
Dental X-rays 
MRI 

fMRI, Fundus camera 
DSA 

Endoscopy, fundus camera 
EM microscopy 
Surgery, radiotherapy 

CT/PET/MRI 



All 



MRI, Microscopy (cytology) 
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The B-spline of degree n satisfies the two-scale relation 
cp n (x/2)=VTIh n (k)cp n (x-k) 

k £ z 

where h n (k) is the binomial filter of order n+l whose trans- 
fer function is 

h n (k)^H"(z)=V2(^^j 

In 1946, Schoenberg proved that any polynomial spline 
of degree n with knots at the integers could be represented 
as a linear combination of shifted B-splines. Thus, our ba- 
sic spline space Vo n can also be specified as 

V;=|s 0 (x) = Ic(k)(p n (x-k)|c6l 2 | 

where the weights c(k) are the so-called B-spline coeffi- 
cients of the spline function s 0 (x). In addition, it can be 
shown that the Bsplines {<p n (x-k)} k£z constitute a Riesz ba- 
sis of Vo" in the sense that there exist two constants An >0 
and Bn< + infinity such that 




The lower inequality implies that the B-splines are 
linearly independent (i.e., s o (x)=0— >c(&)=0). The upper 
inequality guarantees that Vo n eL r Hence, any polyno- 
mial spline has a unique representation in terms of its B- 
spline coefficients c (k). Schoenberg also proved that the 
B-splines are the shortest possible spline functions. This, 
together with the fact that these functions have a simple 
analytical form, makes the B-spline representation one of 
the preferred tools for the study and characterization of 
splines. 

Mathematical model of Bi-orthogonal Wavelets 

In the most general case, the construction of biorthogo- 
nal wavelet bases involves two multiresolution analyses of 
L2: one for the analysis, and one for the synthesis 7. These 
are usually denoted by 

m{(cp~)}. cz and m{ (cp)}. cz 
where <p~(x) and <p(x) are the analysis and synthesis scal- 
ing functions, respectively. Note that (p~(x) and <p(x) can 
be arbitrary solutions of a two-scale relation and not nec- 
essarily the causal B-splines <pn defined previously. The 
corresponding analysis and synthesis wavelets \|/~ (x) and 
\\i(x) are then constructed by taking linear combinations of 
these scaling functions. 

\j> (x / 2) =V2T k g( k )^( x " k ) 



y (x/2)=V2 Z k g(k)cp(x-k) 

They form a biorthogonal set in the sense that 

< ^,k^f, 1 > = § 1 -f,k- 1 
with the short form convention V|/ = 2 ,/2 \|/ (2' x-k). T his al- 
lows us to obtain the wavelet expansion of any L2 -function 
as 

VfeL, f=I I (f,¥ 1;k )¥ 1(k 

1 iez kez 
Note that the underlying basis functions are usually 
specified indirectly in terms of the four sequences h(k), 
h~(k), g(k) and g~(k), which are the filters for the fast wave- 
let transform algorithm. 

Mathematical model of Spline Wavelets 

We have a spline wavelet transform whenever the syn- 
thesis functions (\|/(x) and <p(x)) are polynomial splines of 
degree n. This means that the synthesis wavelet can also be 
represented by its B-spline expansion 

y(x/2)=Iw(k)cp°(x-k) 

k€ z 

It is important to observe that the underlying scaling 
function (p(x)eFb n is not necessarily the B-spline of degree 
n — unless h(k) is precisely the binomial filter. This func- 
tion is usually specified indirectly as the solution of the 
two-scale relation 

<p(x/2) =V2"Zh(k)(p(x-k) 

k6Z 

where h(k) is the corresponding (lowpass) reconstruc- 
tion filter. However, in the spline case, there will always 
exist a sequence p(k) such that 

cp(x) = Xp(k)cp n (x-k) 

Such specific B-spline characterizations for various 
kinds of spline scaling functions (orthogonal, dual, or 
interpolating) can be found elsewhere. Note that the se- 
quence p(k) defines an invertible convolution operator 
from 12 into 12 which performs the change from one coor- 
dinate system to the other (i.e., (p to <p«). The basic require- 
ment for tp(x-k) to form a Riesz basis of VnO is that there 
exist two constants Ap>0 and Bp < +cosuch that Ap <| P(e 
'"'/I <B p almost everywhere, where P(e '"' ) denotes the 
Fourier transform of p. 

By combining above <p(x) equation with Bi-orthogonal 
\|/(x/2) equation, we obtain the B-spline coefficients of the 
wavelet \|/(x): 

w (k) = (p * g) (k) W(z) = P(z)G(z) 

Many kinds of spline wavelets have been described in 
the literature. The four primary types can be differenti- 
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ated on the basis of their orthogonality properties; they are 
summarized in Table 3. Four corresponding examples of 
cubic spline wavelets and their duals are also shown in 
Figure 5. 

PROPOSED METHOD 

The method that we proposing in this paper consist of 
three following major steps 

1) De-noise the CT Image using un-decimated Wavelet 
transform; 

2) Compress the CT image using Spline Bi-orthogonal 
Wavelet; 

3) Accurate detection of nodules with the compressed 
image of step 2 using spline wavelets. 

De-noise using un-decimated Wavelet transform 

Developing Image denoising algorithms is a difficult 
task since fine details in a medical image embedding di- 
agnostic information should not be destroyed during noise 
removal. Many of the wavelet based denoising algorithms 
use DWT (Discrete Wavelet Transform) in the decom- 
position stage which is suffering from shift variance. To 
overcome this in this paper we are proposing the denoising 
method which uses Undecimated Wavelet Transform to 
decompose the image and we performed the shrinkage op- 
eration to eliminate the noise from the noisy image. In the 
shrinkage step we used semi-soft and stein thresholding 
operators along with traditional hard and soft thresholding 
operators and verified the suitability of different wavelet 
families for the denoising of medical images. 

Naga Prudhvi Raj and venkateswaralu (7) applied this 



un-decimated wavelet for denoising the leg knee CT imag- 
es and their statistical results are tabulated below in Table 
3. Table has the Mean square error (MSE) and Peak Signal 
to Noise Ration (PSNR) values for different sigma values. 

Hard, soft, semi-soft and stein are different Threshold- 
ing criteria. 

Based on the analysis of the above result we proposed 
to use this Undecimated wavelet approach for de-noising 
the lung CT image before applying the nodule detecting 
algorithm. 

Compression using Spline Bi-orthogonal Wavelet 

With extensive digitization of data and increasing use 
of Computer Aided Diagnosis, most of image data in hos- 
pitals are stored in digital form using picture archiving 
and communication systems. The need for data storage 
and bandwidth requirements is increasing and lossy com- 
pression techniques have become a necessity. The suc- 
cessful use of the wavelet transform in the field of image 
compression has been extensively studied in literature. 
We proposed to use the Biorthogonal spline wavelet (8) 
as an effective mode for medical image compression. The 
research made by Jane P and his team at Thoracic division, 
Department of Radiology, New York University (9) medi- 
cal center achieved the 1:10 ratio of compression without 
any sacrifice in the nodule detection. 

Their results are tabulated as follow: 

The above results make sure that Bi-orthogonal wave- 
let method is suitable for compressing the lung images 
before applying the nodules segmentation technique. In 
the above research they used some filtering methods for 
nodules detection, but in this proposed method we are go- 



Table 3. Classification of Spline wavelets with its main properties 



Wavelet type Orthogonality Compact support Key properties Implementation 



Orthogonal splines 
(Battle-Lemarie, Mallat) 


Yes 


No 


* Symmetry & regularity 
+ Orthogonality 


IIR/IIR 


Semi-orthogonal splines 
(B-splines) (Chui-Wang, 
Unser-Aldroubi) 


Inter-scale 


Analysis or Synthesis 


* Symmetry & regularity 

+ Optimal time-frequency local 
ization 


Recursive IIR/FIR 


Shift-orthogonal splines 
(Unser-Thevenzaz-Aldroubi) 


Intra-Scale 


No 


* Symmetry & regularity 
+ Quasi-orthogonality 
+ Fast decaying wavelet 


IIR/IIR 


Biorthogonal splines 
(Cohen-Daubechies-Feauveau) 


No 


Yes 


* Symmetry & regularity 
+ Compact support 


FIR/FIR 
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SYNTHESIS WAVELET ANALYSIS WAVELET iftx) 




(b) Semi-orthogonal B-spline 




(c) Shift-orthogonal spline (hybrid degree n = 3 and fl = 1) 




t 



(d) Biorthogonal spline (1 = 4, 1=6) 

ure 5. Examples of four different types of cubic splines wavelets and their corresponding duals. 
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ing to use the Spline techniques to detect the nodules in 
the lungs. As we discussed in section 3, the spline wavelet 
techniques are optimum for medical image segmentation 
and object detection. The globe round shaped nodules can 
be efficiently detected using B-spline wavelet. The first 
step of our proposed method Denoising lead to accurate 
detection of lung nodules at step 3 and also the compres- 
sion of images that we proposed in step 2, reduce the com- 
puting time of step 3. Hence in terms of accuracy and time 
this proposed method is optimum. 

CONCLUSION 

Many Nodule detection algorithms were surveyed 
and discussed in this paper. After analyzing the different 
methods, we proposed the optimum method for detecting 
the nodules in lungs which comprise of 3 steps say de- 
noising, compression and Spline based nodule detection 
respectively. It has been proposed to use Un-decimated 
wavelet for Denoising, Bi-orthogonal wavelet method for 
Compression and B-Spline Wavelet technique for Nodule 
detection after analyzing many previous research and lit- 
eratures. 
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